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We demonstrate time resolved driving of two-photon blue sideband transitions between super- 
conducting qubits and a transmission line resonator. Using the sidebands, we implement a pulse 
sequence that first entangles one qubit with the resonator, and subsequently distributes the en- 
tanglement between two qubits. We show generation of 75% fidelity Bell states by this method. 
The full density matrix of the two qubit system is extracted using joint measurement and quantum 
state tomography, and shows close agreement with numerical simulation. The scheme is potentially 
extendable to a scalable universal gate for quantum computation. 



In the pursuit of a scalable architecture for quantum 
information processing, Josephson-j unction based super- 
conducting circuits are currently a leading candidate [l[ . 
Controlled coherent two-qubit interactions, essential for 
the realization of universal quantum computation, have 
been demonstrated using a variety of coupling schemes 
@- IE El IE IE S HI • A promising recent advance in the 
field has been the realization of circuit quantum eletro- 
dynamics (QED) [9] in which superconducting qubits are 
strongly coupled to single photons in a transmission line 
resonator. This architecture has been used to couple 
two qubits over a distance of several millimeters by ex- 
change of a virtual photon between two resonant trans- 
ition qubits detuned from the resonator [8[ and by ex- 
change of a real photon between two phase qubits in res- 
onance with the resonator [7] . The possibility to achieve 
such non-local couplings is advantageous for the realiza- 
tion of multi-qubit systems. It is also desirable to be able 
to switch two-qubit couplings on and off with high accu- 
racy, such that high fidelity operations can be carried 
out both on individual qubits, and on specifically cho- 
sen pairs of qubits, avoiding spurious couplings to other 
qubits in the system. This issue has been addressed in 
superconducting circuits by controlling the interaction 
strength between two qubits via an external parameter 
0, IE 13 , or by tuning transitions into and out of reso- 
nance [E 0, m . 

In this Letter we present an experimental realization of 
a two-qubit coupling scheme in circuit QED that makes 
use of sideband transitions between the joint qubit and 
resonator states. In this scheme, qubit transition fre- 
quencies may be kept at values chosen for optimal coher- 
ence, and do not need to be frequency tuned. Since the 
coupling is only present when strong microwave drives 
are turned on to induce the sideband transitions the tun- 
ability or on/off ratio of the coupling can be large. The 
scheme is hence potentially scalable in a similar manner 
to that achieved with sidebands in trapped ion systems 

Til, El. 



In the dispersive regime of circuit QED, the qubit tran- 
sition frequency c^a is detuned from the resonator fre- 



quency oj t by A = \uo v — uo a\ ^> 9, where g is the cou- 
pling strength between resonator and qubit on resonance 
(see Fig. Ufa)). Although the qubit and resonator do 
not directly exchange energy in this case, the residual 
dispersive coupling still allows sideband transitions link- 
ing the qubit and resonator states to be accessed using 
strong additional microwave drive fields [IE Hi]. Blue 



sideband transitions involve the simultaneous excitation 
of both qubit and resonator, at a transition frequency 
u\ = u r + cja, while the red sideband involves the ex- 
change of an excitation between the two systems, at a 
transition frequency uo^ = \uj t — ua\- Such sideband tran- 
sitions were first observed in a superconducting circuit in 
an experiment coupling a flux qubit to a superconducting 
quantum interference device (SQUID) [l5j |. 

Due to symmetry considerations in the circuit QED 
Hamiltonian, single photon sideband transitions with ei- 
ther a Cooper Pair Box (CPB) [IE biased at charge de- 
generacy, or with a transmon qubit [IE EE , are forbid- 
den to first order. They may however be accessed us- 
ing two photons, whose sum or difference frequency is 
equal to one of the sideband transition frequencies uj^J 
[141 ]. A spectroscopic investigation of these two-photon 
transitions for a CPB has been reported previously, in 
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FIG. 1: (a) Schematic of the energy levels of a harmonic oscil- 
lator with resonant frequency u r , and levels \n) corresponding 
to photon number n, and qubit with transition frequency cja, 
and ground and first excited states \g) and |e) respectively, (b) 
Combined level diagram, indicating the blue sideband transi- 
tion uo~^ = oo r + cja, which is forbidden to first order, but may 
be driven using two photons, for example both at a frequency 
of cj+/2. 
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FIG. 2: (a) False color optical microscope image of the two- 
qubit circuit QED device. The sapphire substrate is shown 
in dark green, the Niobium resonator in light blue, and the 
Aluminium qubits in dark blue, (b) High magnification false 
color optical microscope image of one of the qubits coupled to 
the transmission line resonator, and a local microwave drive 
line, (c) Electrical circuit representation of the two-qubit cir- 
cuit QED device with local controls. The magnetic fluxes <J>a 
and <E>b through the SQUID loops of the qubits are tuned us- 
ing local superconducting coils (green), while transitions are 
driven using local microwave drive lines (blue). 



which two microwave drive fields of different frequency 
were applied to the resonator, detuned slightly from the 
qubit and resonator frequencies [l3j]. The increased cou- 
pling strength of the transmon to the resonator [17] when 
compared to the CPB allows driving of the sidebands at 
larger detunings of the drives from the qubit and res- 
onator. Here we choose to drive the blue sideband with 
two photons of equal energy, using a single microwave 
drive of frequency c^^/2 = (uo T + ooa)/2 (see Fig. OJb)). 
In this configuration the drive is equally detuned from 
the qubit and resonator, maximizing the selectivity of 
the sideband transition with respect to undesired off- 
resonant driving of the bare qubit transition and off- 
resonant population of the resonator. 

We work with a two qubit circuit QED device, for 
which an optical image and circuit schematic are depicted 
in Fig. [21 A transmon qubit is fabricated at each end of a 
coplanar waveguide resonator with bare fundamental res- 
onant frequency cj r /27r = 6.44 GHz. The qubits, labelled 
hereafter as qubit A and B are near identical, with a cou- 
pling strength to the resonator g/2n = 133 MHz and a 



charging energy Ec/h = 232 MHz. The transition fre- 
quencies of the qubits are tuned independently by chang- 
ing the local magnetic fluxes <£a and $b through the 
SQUID loops of qubits A and B respectively, using two 
small superconducting coils mounted below the chip. Di- 
rect qubit and sideband transitions are driven selectively 
using microwave transmission lines coupled capacitively 
to each qubit locally (see Fig. Efc)). For the experi- 
ment, the qubits are tuned to transition frequencies of 
loa/^tt = 4.50 GHz and cjb / 27r = 4.85 GHz respectively, 
well into the dispersive regime. 

Sweeping the frequency of a single strong microwave 
drive on one of the selective qubit drive lines, and si- 
multaneously measuring the resonator transmission to 
carry out a dispersive measurement of the qubit state 
[l9t shows both the direct qubit transition at uja, and an 
additional spectral line corresponding to the two-photon 
blue sideband at uo\/2 (see Fig. [3{a)). The high selec- 
tivity of the individual qubit drive lines suppresses the 
coupling to qubit B, the transitions of which are barely 
visible in Fig. [3ja). The sideband spectral line is shifted 
down at high drive amplitudes due to the AC- Stark ef- 
fect on the bare qubit transition [20j|. At sufficient drive 
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FIG. 3: (a) Spectroscopy on the selective drive line of qubit 
A, showing the direct transition of the qubit at 4.50 GHz, 
and the two-photon blue sideband at 5.47 GHz. (b) Time re- 
solved measurements of the excited state population of qubit 
A as a function of the length of a blue sideband pulse. Blue 
solid circles show the response for the initial state |#0), and 
green open circles for |e0). Solid lines: Master equation simu- 
lations using decay parameters extracted from separate mea- 
surements. (Inset) Level diagram showing competing drive 
and decay rates in the system. 
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amplitude, time resolved Rabi oscillations are observed 
on the blue sideband transitions (see Fig. 09(b)) by ap- 
plying microwave pulses of fixed amplitude and varying 
length, at the Stark-shifted two photon blue sideband 
frequency, to the local drive line of one qubit. The qubit 
excited state population is then measured by applying 
a microwave pulse to the resonator, and comparing the 
time evolution of the transmission to a numerical solution 



of the coupled qubit resonator dynamics [21]. 

Fig. [3](b) shows the extracted excited state population 
for such an experiment on qubit A (average of 6.6 x 10 4 
repetitions). The solid blue data points show the re- 
sults obtained when the system is initially in its ground 
state |g0), while the green open data points correspond 
to an experiment in which the system is first excited to 
the state |e0) using a 7r-pulse on the direct qubit transi- 
tion. In the latter case, no oscillations are observed since 
the state |e0) is not addressed by the blue sideband, a 
fact that was also observed in Ref. [TH. Also shown in 
Fig 09(b) (solid lines) are master equation simulations of 
the evolution, using values for the qubit relaxation rate 
7/27T = 0.2 MHz and photon decay rate k/2tt = 1.7 MHz 
taken from separate measurements, and with the blue 
sideband transition rate Q/2tt = 8.15 MHz as a fit pa- 
rameter. 

At long times it can be seen in Fig. 09(b) that the 
qubit excited state population tends to a steady state 
value of p ss « 0.9. This can be explained by con- 
sidering the different drive and decay channels present 
in the system (see inset). Since the state |e0) has no 
blue sideband transition, and 7 <C ft, driving the side- 
band from the ground state \g0) results in a build up 
of population in |e0). A simple model taking into ac- 
count only the lowest four levels of the system predicts 
p ss — + 7ft + 7^) = 0.87. It is worth noting that 

an increase in the ratio ft/7 would allow high fidelity ex- 
cited state preparation by pumping of the blue sideband, 
in a similar manner for example to spin preparation by 
optical pumping in semiconductor quantum dots 22|, [23[ . 

Having characterized the sideband Rabi oscillations for 
both qubits A and B, we implement a sequence of pulses 
on the two-qubit plus resonator system to generate qubit- 
qubit entangled states, in a similar manner to that imple- 
mented in trapped ions (24|. The full sequence is shown 
schematically in Fig. [H At time to, the system is in its 
ground state \ggO), and a resonant 7r-pulse is first applied 
to qubit B, generating the state \geO) at t\. A 7r/2-pulse 
on the blue sideband of qubit A then generates the en- 
tangled state (\geO)-\-e 1 ^ \eel))/y/2. The qubit-resonator 
entanglement is then transferred to qubit-qubit entan- 
glement with a 7r-pulse on the blue sideband of qubit B, 
generating the Bell state |^) = (\ge) + e i(f) \eg))/\/2 at t 3 , 
with the resonator returning to its ground state |0). The 
phase (j) of the state is dependent on the phase difference 
between the two blue sideband pulses in the sequence. An 
additional resonant 7r-pulse may now be applied to qubit 
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FIG. 4: Pulse sequence implemented on the two-qubit and 
resonator system, to generate and characterize qubit-qubit 
Bell states. Direct resonant qubit pulses are shown in white, 
while blue sideband pulses between resonator and qubits A 
and B are shown in light and dark blue respectively. The 
ideal state of the system neglecting decay and imperfections 
after step i, at each of the intermediate times t% is indicated 
in the text. 
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A to generate the Bell state |3>) 
at £4. 

A full tomographic reconstruction of the two-qubit 
state is carried out by repeating the pulse sequence 16 
times, with different combinations of additional final sin- 
gle qubit rotations (identity, tt/2 x , 7r/2 y , and ir x on each 
qubit). A microwave pulse close to the resonator fre- 
quency is then used for an averaged measurement of the 
two-qubit state. Since the resonator frequency is shifted 
by a different amount for each of the four two-qubit states 
\gg), \ge} : \eg), \ee) [8], and since the resonator transmis- 
sion quadrature amplitudes depend non-linearly on this 
shift, the full density matrix can be parametrised in terms 
of the 16 averaged measurement results [25|. The closest 
physical density matrices to those experimentally mea- 
sured are found using a maximum likelihood approach 

Examples of experimental two-qubit density matrices 
p measured after carrying out the pulse sequence for 
generation of the = (\ge) + \eg))/y/2 and |$ + ) = 

(\gg) + \ee))/y/2 Bell states are shown in Fig. [51(a) and 
(b). In each case 6.6 x 10 5 repetitions of the experi- 
ment were averaged. The fidelities T = ({^p^)) 1 / 2 of 
the measured states p with respect to the perfect Bell 
states i(j = and $ + are T = 74% and 76%, respec- 
tively, where T is limited primarily by the resonator pho- 
ton lifetime, 1/k « 100 ns. We have repeated the state 
preparation 33 times and calculated the fidelities sepa- 
rately each time, obtaining an average fidelity of 75 ± 2% 
for the full set. The entanglement of formation [27[ of the 
generated states extracted from their density matrices is 
0.09 ±0.04. 

Also shown in Fig. [5] ((c) and (d)) are the density ma- 
trices resulting from a master equation simulation of the 
experiment. The simulation makes both the dispersive 
and rotating wave approximations. Two energy levels 
for each qubit, and the first 5 levels of the harmonic os- 
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FIG. 5: (a,b) Real and imaginary parts of the two-qubit density matrix extracted from two-qubit tomography measurements 
after carrying out the pulse sequence to generate the Bell states |^+) and |3>+) respectively. (c,d) Simulated density matrices 
corresponding to the states shown in (a) and (b) respectively. 



cillator are taken into account. Energy relaxation and 
decoherence rates for the three individual systems are 
taken from separate measurements. The fidelity between 
simulated and measured density matrices is 97 ± 2%, in- 
dicating that the loss of fidelity in the experimental state 
preparation is dominated by photon loss. 

The excellent agreement between experiment and sim- 
ulation gives confidence that with a longer photon life- 
time, higher fidelity entangled states are likely to be 
within reach by this method. The photon lifetime may 
be increased in future experiments with the use of higher 
quality factor resonators. The scheme may then be ex- 
tendable to the implementation of a universal CNOT gate 
[ill . [IK , and should scale well to a system with a larger 
number of qubits. The sidebands may also be used to 
generate non-classical states of the microwave field in the 
resonator, such as Fock states. 
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